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Abstract 

We investigate two-loop gluino corrections to the effective Lagrangian for b — > s + 7(5) in 
the minimal supersymmetric extension of the standard model (MSSM) at large tan /3, including 
the contributions in which quark flavor change is mediated by charginos. Using the translation 
invariant of loop momenta and the Ward-Takahashi identities (WTIs) that are required by the 
SU(3) C x U (l) em gauge invariance, we simplify our expressions to concise forms. As an example, 
we discuss two- loop gluino corrections to the CP asymmetry of inclusive B — > X s "f decay in CP 
violating MSSM. 
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I. INTRODUCTION 

The measurements of the branching ratios at CLEO, ALEPH and BELLE 1] give the 
combined result 

BR(B -> Ij) = (3.11 ± 0.42 ± 0.21) x 10" 4 , (1) 
which agrees with the next-to-leading order (NLO) standard model (SM) prediction 

BR(B -f X s7 ) = (3.29 ± 0.33) x 10^ 4 . (2) 

Good agreement between the experiment and the theoretical prediction of the SM implies 
that the new physics scale should lie well above the electroweak (EW) scale. The systematic 
analysis of new physics corrections to B — > X s 7 up to two-loop order can help us under- 
standing where the new physics scale sets in, and the distribution of new physical particle 
masses around this scale. In principle, the two-loop corrections can be large when some ad- 
ditional parameters are involved at this perturbation order beside the parameters appearing 
in one loop results. In other words, including the two-loop contributions one can obtain a 
more exact constraint on the new physics parameter space from the present experimental 
results. 

Beside the Cabibbo-Kobayashi-Maskawa (CKM) mechanism, the soft breaking terms 
provide a new source of CP and flavor violation in the MSSM. Those CP violating phases 
can affect the important observables in the mixing of Higgs bosons [sj, the lepton and 
neutron's electric dipole moments (EDMs) jj, Q], lepton polarization asymmetries in the 
semi-leptonic decays jfj], the production of P-wave charmonium and bottomonium and 
CP violation in rare 5-decays and in B°B° mixing Q|. At present, the strictest constraints 
on those CP violation phases originate from the lepton and neutron's EDMs. Nevertheless, 



if we invoke a cancellation mechanism among different supersymmetric contributions |4j , or 



Q, 

choose the sfermions of the first generation heavy enough j5j, the loop inducing lepton and 
neutron's EDMs bound the argument of the /1 parameter to be < 7r/(5tan/?), leaving no 
constraints on the other explicitly CP violating phases. 

The supersymmetry models at large tan f3 are implied by grand unified theories, where the 
unification of up- and down- type quark Yukawa couplings is made From the technical 
viewpoint, the dominant contributions to the relevant effective Lagrangian are the terms 
proportional to (tan/3) n (n =1, 2, • • •) in a large tan/3 scenario. This will simplify our 
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two-loop analysis drastically since we just keep those terms enhanced by tan (3. 

Assuming no additional sources of flavor violation other than the CKM matrix elements, 
the authors of Q present an exact ana, y sis of the two,oop g luino c—as to the rare 
decay b — > s + 7(g) in which quark flavor change is mediated by the charged Higgs in CP 
conserving MSSM at large tan (3. They also compare their exact result with that originating 
from the heavy mass expansion (HME) approximation Although the HME result 

approximates the exact two-loop analysis adequately when the supersymmetry energy scale 
is high enough, their analysis implies that the difference between the HME approximation 
and exact calculation is obvious in some parameter space of the MSSM. However, they 
do not consider the case in which quark flavor change is mediated by the charginos (the 
super partners of the charged Higgs and W bosons). In fact, we cannot provide any strong 
reason to ignore the contribution from the diagrams in which quark flavor change is induced 
by the charginos, even within large tan/5 scenarios. In this work, we present a complete 
analysis on the two-loop gluino corrections to the rare transitions b s + r y(g) by including 
the contributions of those diagrams where quark flavor change is mediated by charginos 
in the framework of CP violating MSSM at large tan/3. Furthermore, we also simplify 
our expressions to concise forms through loop momentum invariant and the WTIs that are 
required by the SU(3) C x U(l) em gauge invariance. 

The paper is organized as follows. In Sec. |Hj we give all the diagrams needed to evaluate 
the 0(a s tan (3) contributions to the Wilson coefficients C 7 and C 8 entering the branching 
ratio BR(B — > X^). The corresponding Wilson coefficients at the matching EW scale /i EW 
are also presented there. We apply the effective Lagrangian to the rare decay B — ► X s 7 in 
Sec. IIHI By the numerical method, we show the two-loop corrections on the CP asymmetry 
for the process. Our conclusion is given in Sec. IIV[ and some long formulae are collected in 
appendices. 

II. THE WILSON COEFFICIENTS FROM THE TWO-LOOP DIAGRAMS 

In this section, we derive the relevant Wilson coefficients for the partonic decay b — *■ sj 
including two-loop gluino corrections. In a conventional form, the effective Hamilton is 
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(a) (b) 

FIG. 1: The self energy diagrams which lead to the magnetic and chromo-magnetic operators in 
the MSSM, the corresponding triangle diagrams are obtained by attaching a photon or gluon in 
all possible ways. 



written as 

H n = -^v;x b tc^)o i} (3) 

where V is the CKM matrix. The definitions of the magnetic and chromo-magnetic dipole 
operators are 

= j^y 2 m b ^)s L T a ^ v b R Gl , (4) 

where and G a v are the field strengths of the photon and gluon respectively, and T a (a = 
1, 8) are SU(3) C generators. In addition, e and g g represent the EW and strong 

couplings respectively. The other operators O i (i — 1, ■ • • , 6) are defined in 

In the framework of CP violation MSSM, the one-loop analysis on the CP asymmetry of 
inclusive B — > X s 7 decay has been presented elsewhere The two loop gluino diagrams, 
contributing at (9(a s tan/3) to the Wilson coefficients of the magnetic and chromo-magnetic 
dipole operators, are obtained from the self energy diagrams (a), (b) of FIG. ^ by attaching 
a photon or gluon in all possible ways. The calculation of the Wilson coefficients for the 
operators in Eq. (j3J) at the two loop order is more challenging than that at the one loop 
order. Before we give those Wilson coefficient expressions explicitly, we state firstly the 
concrete steps required to obtain the coefficients from those two loop diagrams. 

• After writing the amplitudes of those two loop triangle diagrams, we expand them in 
powers of external momenta to the second order. 

• The even rank tensors in the loop momenta qi, q2 can be replaced as follows 
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d D q 1 d D q 2 qi^Qiv, qi^u 9^ f d D q 1 d D q 2 q\, q x ■ q 2 



(2n) D {2ti) d V D J (2n) D (2tc) d V 

d D q 1 d D q 2 qi^qiuqipqia, q\^q\vq\ P q2u 
(2tt) d (2tt) d V 

T u vpa f d D( li dD( l2 qf, qf(qi • 92) 



D{D + 2)J {27i) d {2ti) d V 
d D q 1 d D q 2 q^q^pq^ 
(2tt) d (2tt) d V 

I d D q 1 d D q 2 1 , D{qy ■ q 2 ) 2 -gfgf (gi • g 2 ) 2 - gigg \ 

~~ y (27r) D (27r) D D q V £>(£>- 1)(D + 2) M,ypCT D(D-l) 9 ^ 9p ^ ' U 

where D is the time-space dimension, T ^ ct = + g„ p g„„+ g^g pv and Z> = ((g 2 - 

qi) 2 —rn 2 )(q 2 —m 2 )(q2—m 2 ). The odd rank tensors in the loop momenta can be dropped 
since the integrations are symmetric under the transformation q\p — > —q\ i2 . Here, we 
only retain the simplest two-loop propagator composition 1/T> which corresponds to 
the two-loop vacuum diagram (FIG. |2J). Any complicated composition of two-loop 
propagators can be expressed as the linear combination of the simplest one 1/T> by 
use of the obvious decomposition formula 

1 1 / 1 1 N 

(Q 2 — m 2 )(Q 2 — m 2 ) m 2 — m 2 B ^Q 2 — m 2 Q 2 — m 2 B > 
with Q = qi, q 2 , or q 2 — q 1 . As an example, we apply the above two steps to the triangle 
diagram in which a external photon is attached to the internal squark s i line (FIG. 
El). After expanding the corresponding amplitude in powers of external momenta to 
the second order, we have 
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FIG. 3: A triangle diagram in which the external photon is attached to squark s v 



iAl(p,k) = ^-e^^^K^^tan/?)^)^^)^^^^^-^^) 



w w 

v / 2m w s w A r rf^ d D g 2 1 

e 1 {2n) D {2Tt) D V H {ql-m 2 ) 



1 q{ — m. q 2 — m z > 



H+ 



x{(Zt). 3 i4 2 {2q 1 -2p-k)^ + 

-m t \m 3 \e w * (zj) . a (2 9l - 2p - A;) M cu + } , (7) 

where p, fc represent the incoming momenta of the external quark b and photon 
respectively, V H = ((q 2 - qif -m 2 )(q{ - \m 3 \ 2 )(ql -m?.)(gf -m 2 ){q\ -m 2 H+ ) (h 3 = 
1, 2) and e d = —1/3. Z~ (q = u, d, • • • , t) are the mixing matrices of scalar quarks, 
and A q are the corresponding trilinear soft breaking parameters. Furthermore, 9 3 
denote the CP phases of the SU(3) C gaugino mass and of the fx parameter respectively. 
Since the quark mass m b from the Yukawa coupling of bottom quark is same order as 
the external momenta p, k in magnitude, we just expand the propagators in powers of 
the external momenta to the first order. In the soft breaking potential, the CP phase 
9 3 is contained in the gluino mass terms 

K|e^A G A G + K|e-^A G A G , (8) 
where A G denotes the gluino in a two-component Majorana spinor. With the redefini- 
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tion of the gluino field 



(9) 



the mass terms are transformed into 



\ m 3 \99 



(10) 



with the four-component Majorana spinor 
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(11) 



Correspondingly, the CP phase 8 3 is transfered from the mass terms to the quark- 
squark-gluino vertex which is given by [13] 



H.c. 



(12) 



Here, a, [3 = 1, 2, 3 are quark and squark color indices, and u± = ^ 2§ -- This is the 
reason why there is a |m 3 | rather than m 3 in the gluino propagator. Using Eq. Eland 
Eq. El the amplitude of FIG. El is finally formulated as 



iAl(p,k) = -i- e ^_^(^-tan ) 9)(Z I ) ai4 (Z;) <ia (| / iKe^(Z f ; 



2,J 



d D qi d D g 2 



4 g^i • g 2 



D gf — m 2 
i ,D{q x -q 2 y -q\q\ 



(2p + k)fj,u+ -qi-q 2 (2p + k) f 



ql-ml + \ D(D-l) 



P»» + ~ D{D _ l) 7>H 



-m t |m 3 |e i9 3(zt) 



J',2 



D q( — m? 



D q2 — m 



(2p + fc)„w+] } . 



H+ 



(13) 



In a similar way, we can obtain the other triangle diagram amplitudes. 
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Using loop momentum translation invariant, we formulate the sum of those amplitudes 
in gauge invariance form explicitly, then extract the corresponding Wilson coefficients 



which are expressed by the two- loop vacuum integrals [14|. In fact, there are many 
identities among those two-loop integrations. In order to obtain those necessary iden- 
tities which are used to simplify the sum of those triangle amplitudes, we start from 
the zero integration such as 

d D q 1 d D q 2 q\ -p i\{ 



Under the loop momentum translation q 2 — > q 2 —a where a is same order as the external 
momentum p in magnitude, we expand the integration in powers of the momenta p, a 
to the second order 

dPq x d D q 2 q x ■ p , , , 



x 



(2n) D (2ir) D V 
d D q 1 dPq 2 q\ ■ p 
(2tt) d {2tt) d V 

( h - ii) - ii 

d D qi d D q 2 1 



2(?2 - qi) 



(?2 - 1l) 



+ 



2q 2 ■ a 



(2tt) d (2tt) d V 
qf - 2qjq x ■ q 2 



D 



)(p-a) 



q\,, , 2 
D W + (q 2 - qi ) 2 - ml 
(?i • <?2) 2 - q\q\ 



D{qi-q 2 f 



q\q\ 



D(D-l) 



D(D-l) 



2 U D(qi-q2) 2 

D(D - 

(qi ■ q 2 ) 2 ~ q\q\ 



q 2 — m 2 



q\q\ 



V 



gigi • q2 

D 



){p-a) 



D(D-l) 
The above identical equation implies 
d D qi d D q 2 1 r 1 







(15) 



(2n) D (2ir) D V I (q 2 - qi f - m 2 1 D(D — 1 
1 ,D(q 1 ■ q 2 f - q\q\ q\q x ■ q 2 



D(qi ■ q2? ~ q\q\ qf - 1q\q\ • <h y 



D 



m? a 1 D(D 
d D qi d D q 2 1 



1: 



D 



= 



(qi ■ q2? - q\q\ 



(27r) D (27r) D V t(q 2 - qi ) 2 -mt 

2 (gi ■ g 2 ) 2 - q\ql , 9i 1 = 
g|-m 2 D(D-l) £>J" 



. 



(16) 



Similarly, we can get the following identities from the invariant of Eq. El under the 
loop momentum translation q% — > gi — a, q 2 ^ q 2 — a: 



d D qi d D q 2 1 r 2 + L> 



(27r) B (2tt) d X? c 



{ 



-gi ■ (?2 



gi + — 
qf 



2 gfgi • (52 - Qi) 



D 
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q\ — ^ 2 



D(qi ■ g 2 ) 2 - q\q\ q\q\-qr , _ () 



D(D-l) D 

dPqx d D q 2 1 r gi • (g 2 - gi) _ 2 fa ■ g 2 ) 2 - g 2 g 2 . = 

(2tt) d (2tt) s I? 1 D g 2 -m 2 D(£> - 1) J 1 J 



Using the concrete expressions of two-loop vacuum integrals in Ref. 14], we can also 
verify those equations in Eq. E3and Eq. El directly after some tedious calculations. 
Replacing the numerator of Eq. El with other odd rank tensors in the loop momenta 
qi, q2, the more additional identities among two- loop integrations are gotten. In gen- 
eral, those identities are linearly dependent. After some simplification, we obtain those 
linearly independent equations in appendix [X] Certainly, those linearly independent 
equations can also be derived from those two-loop integrations in which the numera- 
tors are even rank tensors of the loop momenta qi, q%. However, the process to derive 
the linearly independent equations with the numerators in even powers of the loop 
momenta is more complicated than that with the numerators in odd powers of the 
loop momenta. 

After the above procedure, we finally obtain the relevant coefficients from the charged 
Higgs contribution up to 0(a s tan (3) 



CU/O = ^(47r) s ei (a.tan^(^ w (2j) 4ia (|/*Ke- <tf M(Z I ) !U 



(2tt) 4 (2tt) 4 D h 
-rn t \m 3 \e^{Z%Ml^ , 



3 

V2s w m v A \ f _^3L_^lJ_( (z U u-g 

f e WW,/ (27r) 4 (27r) 4 P H t ( ^^-^hw 

^ 3 (^), 2 A^ {2) }, (18) 



— mjm, 



where a g = q 2 /47r. and the expressions of the form factors J\f 1, 9 can be found in appendix 
IB! Note, Ref. 10] has also obtained the Wilson coefficients from the same diagrams. We 
formulate our expressions in the more concise forms using the identities from appendix [XJ 
For the chargino contribution that is ignored by Ref. [10j, we can similarly have 
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Cr^OO = ^e d (4 7 r) 3 ( as tan/5)(Z.) 2 , i (Zt) i2 (Zt) jl (Zi) fei: 
+v / 2m w m t (Z.) 2 ,(Z_) 1 ,^ 

X fc (4) 

C 8 ,M = ^(47r) 3 ktan/3)(Z,) 2 , i (2t) i]2 (Zt), 1 (zl) M 



3 

y Ai A2 1 r , . , , rg 

7 (2vr) 4 (2vr) 4 D Xfe l * XfcV ^ 2jV +;2 ' fc x±d) 

+v / 2m w m t (Z.) 2 ,(Z_) 1 ^ 

+m>J m3 | e ^(2 { ) y (2 + ) 2 ^ } , (19) 

with r> Xfc = ((92 — <?i) 2 — rf)(<li ~ \ m -j,\ 2 ){(li — m 1 )(<?2 — m ~.)(?2 — m x fc )' ^-> ^+ are the left- 
and right-handed mixing matrices of charginos, and the form factors A/" 7 ; 9 (i — 1, 2, 3, 4) 
are collected in appendix El After we simplify the sum of the sby (g) triangle diagram 
amplitudes using the identities in appendix |X1 we find that the effective sbj (g) vertices 
should also include the two-point operator 

o sc = j^ 2 mMs L m 2 K , (20) 

beside the magnetic (chromo-magnetic) dipole operators 7 {O a ). Here, the covariant 
derivative acting on the quark fields is 

D^d^-iee^-ig.G^, (21) 

with G M = G a T a . Certainly, the Wilson coefficient of this operator does not give any 
contribution to the rare process b — > 57 after we evolve the corresponding coefficients from 
the matching EW scale to the hadronic scale. Nevertheless, when we extract the Wilson 
coefficients of 7 (O a ), it makes sense to keep this operator for the following reason. Beside 
the effective vertex with two quarks 
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the operator O sc can also induce the effective vertices with two quarks and one photon or 
gluon 

°" ~ (4^ Tam >( 2 ^ + ^ 7 "H (23) 
in the momentum space. Here, p, k are the incoming momenta of the external quark b 
and gauge boson (7 or g) respectively. For the effective shy (g) vertices A 1 ^ (p, k) , the 
corresponding WTIs required by the SU(3) C x U(l) em gauge invariance are written as 

iee d (p(p + k) — £(p)) = ik ■ A y (p, k) , 

igT a (p(p + k) — E(p)) = zA; • A*(p, k) , (24) 

where ££(p) represents the sum of amplitudes for the self energy diagrams (FIG. Q]). Ex- 
panding the self energy amplitudes in powers of external momentum to the third order, we 
have 

= J^y B o m b P 2 ' ( 25 ) 

where the function B only depends on the heavy freedoms which are integrated out. In the 
effective theory, there are two deductions from the WTIs: 

• the effective sb'j (g) vertices can be formulated as 

k) = j£y 2 ee d m b {B^2p^ + frj + B][ft, % ]}u + , 

iAl(p,k) = ^gT a m b {Bf(2p^ + ^)+B^, % ]}u; + (26) 

after we expand ^ (p, k) in powers of the external momenta to the second order, 
where BJ'® are the functions of heavy freedoms only; 

• additionally, B = BJ = B{. 

The above deductions can be taken as the criterion to test our calculations. In Eq. 
the functions B^ 9 are proportional to the Wilson coefficients of the magnetic and chromo- 
magnetic dipole operators respectively. 

As an application, we will investigate the CP asymmetry of the rare decay B — > 
within the framework of MSSM. 
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m„ t (GeV) 

FIG. 4: The CP asymmetry and branching ratio of the inclusive B — > X s j decay versus the 
charged Higgs mass fn H+ ■ Dash-line: theoretical prediction at the one-loop order, and solid-line: 
theoretical prediction at the two-loop order, when tan/? = 30, |m 2 | = 300 GeV, rn 3R = 500 GeV, 
the other parameters are taken as in the text. The gray band is the experimental allowed region 
for the branching ratio BR{B — > X g j) at la deviation. 



III. DIRECT CP VIOLATION IN B -> V s7 



In the SM, the CP asymmetry of the B — > X s 7 process 

A IB : X -A - TiB - - Y - 7) " TiB - X - l) (97) 
A CP (B - X.y) - r(B ^ A . i7)+r(B ^ Xi7) (27) 

is calculated to be rather small: A CP ~ 0.5% [15]. For experimental data, the recent 
measurement of the CP asymmetry implies the 95% range of 

-0.30 < A CP (B -> X s7 ) < 0.14 . (28) 

In other words, studies of the direct CP asymmetry in B — > X s 7 may uncover new sources 



of the CP violation which lie outside t 



prediction has been presented in Ref. 







le SM. Up to the NLO, the complete theoretical 



In order to eliminate the strong dependence on 



the b-quark mass, the branching ratios is usually normalized by the decay rate of the B 
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m„ t (GeV) 

FIG. 5: The CP asymmetry and branching ratio of the inclusive B — > X s j decay versus the 
charged Higgs mass fn H+ ■ Dash-line: theoretical prediction at the one-loop order, and solid-line: 
theoretical prediction at the two-loop order, when tan/? = 60, |m 2 | = 300 GeV, rn 3R = 500 GeV, 
the other parameters are taken as in the text. The gray band is the experimental allowed region 
for the branching ratio BR{B — > X g j) at la deviation. 



meson semileptonic decay: 



F(B -> X 7) _ 6a 
T(B X c ev) ~ ^f(z) 



V*V 



(29) 



where f(z) = 1 — 8z + 8z 3 — z 4 — 12z 2 lnz is the phase-space factor with z = (m c /m b ) 2 , 
and a = e 2 /(4ir) is the electroweak fine-structure constant. The CP asymmetry in the rare 
decay B — > X s 7 is correspondingly formulated as [3| 



a cp (b - x s7 ) = ^#4{^ Im [ c, 2 (^J c ';(^)] - ^Mc.CaOc;^ 



ts t6 



(30) 



with #(2) = (5 + In 2 + In 2 2 - tt 2 /3) + (ln 2 z - 7r 2 /3)^ + (28/9 - 4/31n^)^ 2 + 0(z 3 ). 
The C 2 (/i b ) is the Wilson coefficient of the operator 2 = s L j fl q L q L 'y^b L (q = c, u) at the 
hadronic scale. From now on we shall assume the value BR{B — > X c eu) = 10.5% for the 
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FIG. 6: The CP asymmetry and branching ratio of the inclusive B ^ X s j decay versus the 
parameter \m 2 \. Dash-line: theoretical prediction at the one-loop order, and solid-line: theoretical 
prediction at the two-loop order, when tan/3 = 30, rn R+ = 600 GeV, rn aR = 200 GeV, the other 
parameters are taken as in the text. The gray band is the experimental allowed region for the 
branching ratio BR(B — » X s ^f) at la deviation. 



semileptonic branching ratio, a s (m z ) = 0.118, a(m z ) = 1/127. For the standard particle 
masses, we take m t = 174 GeV, m b = 4.2 GeV, m w = 80.42 GeV, m z = 91.19 GeV and 
z = m c /m b = 0.29. In the CKM matrix, we apply the Wolfenstein parameterization and set 
A = 0.85, A = 0.22, p = 0.22, r] = 0.35 [l^. Without loss of generality, we always assume 
the supersymmetric parameters /x = A t e~ l7T ^ 2 = 100 GeV, m 3 e"'^ 4 = 300 GeV, A 3 e~ l7T ^ 2 = 
m tR = 200 GeV, m tL = m SL = 5 TeV here. In order to suppress the one-loop EDMs, we 
choose the /x parameter CP phase 8 =0. As for the CP phase which is contained in the 
SU(2) gaugino mass parameter m 2 , it is set as 8 2 = arg(rn 2 ) = n/4. 

Taking tan/5 = 30, \m 2 \ = 300 GeV, m SR = 500 GeV, we plot the CP asymmetry and 
branching ratio of the inclusive B — > X s 7 decay versus the charged Higgs mass in FIG. 
@] Considering the experimental constraint on the branching ratio BR{B — > X s ^) at la 
tolerance, the CP asymmetry including the two-loop corrections can be larger than 3%, and 
the one-loop result is smaller than 1% with our chosen parameters. The choice of parameter 
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|mj (GeV) 

FIG. 7: The CP asymmetry and branching ratio of the inclusive B — > X s 7 decay versus the 
parameter |m 2 |. Dash-line: theoretical prediction at the one-loop order, and solid-line: theoretical 
prediction at the two-loop order, when tan/3 = 60, m + = 600 GeV, rn aR = 200 GeV, the other 
parameters are taken as in the text. The gray band is the experimental allowed region for the 
branching ratio BR(B — » X s ^f) at la deviation. 

space of FIG. El is identical with that of FIG. 0] except for tan/3 = 60. After including the 
two-loop corrections, we find that the CP asymmetry can reach 5% with an increasing of 
the charged Higgs mass when tan/3 = 60, while at the same time keeping the branching 
ratio BR(B — ► X^) is within the la deviation experimental bound. Since the two-loop 
correction is proportional to tan /3, we can understand why the differences between the one- 
and two-loop predictions of FIG. El (tan (3 = 60) are larger than that of FIG. 0] (tan j3 = 30). 
From the numerical analysis, we find that the two-loop corrections to the Wilson coefficients 
are still rather smaller than the one-loop results at tan j3 = 30, while the two-loop corrections 
are comparable with the one-loop results at tan/3 = 60. Since the bottom quark Yukawa 
coupling approximates to 1 at tan/3 = 60, it should be argued whether or not we can safely 
apply the perturbative expansion to give the theoretical predictions of physics observables 
for such high tan j3. 

Now, let us study the variance of two-loop results with the soft 577(2) gaugino mass 
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parameter \m 2 \. Taking tan/5 = 30, m SR = 200 GeV, m H+ = 600 GeV, we plot the 
theoretical predictions for the CP asymmetry and branching ratio of the inclusive B — > X s 7 
decay versus the parameter |m 2 | in FIG. [Hjat one- and two-loop order respectively. If the 
theoretical prediction for the branching ratio satisfies with the experimental bound at la 
deviation 

2.48 x 10" 4 < BR(B -> X s7 ) < 3.74 x 10" 4 , 

the CP asymmetry including the two-loop corrections is about ~ 1.5%. The choice of 
the parameter space in FIG. is identical with that of FIG. El except for tan/3 = 60. In 
this scenario, the two-loop prediction on the asymmetry is about ~ 4%. Note that the 
dependence of the two- loop corrections on the parameter |m 2 | is milder than that on the 
charged Higgs mass rn H+ ■ This fact can be understood as follows: the amplitudes of the 
corresponding triangle diagrams depend on the charged Higgs mass in form 1 / (Q 2 — rn 2 ^ + ) 
(Q denotes loop momenta q% , q 2 , or or q 2 — qi), and depend on the parameter \m 2 \ through 
the chargino propagator (Q — m x )/(Q 2 — m 2 ) (m x denotes the chargino mass) before the 
loop momentum integration. 

In our analysis, we do not compare the exact two-loop analysis with the HME result since 
the discussion has already been presented in Ref. jlo| . 

IV. CONCLUSIONS 

In this work, we present the complete two-loop gluino corrections to inclusive B — > 
decay in explicit CP violating MSSM within large tan/3 scenarios. Beside the diagrams 
where quark flavor change is mediated by the charged Higgs, we also include those diagrams 
in which quark flavor change is mediated by the charginos. Using loop momentum transla- 
tion invariant, we formulate our expressions fulfilling the SU(3) C x U(l) em WTIs. From the 
numerical analysis, we show that the two-loop corrections to the branching ratio are com- 
parable with the one-loop predictions at large tan/5. Correspondingly, the CP asymmetry 
can also reach about 5%, which is much larger than that predicted by the SM. 
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APPENDIX A: IDENTITIES AMONG THE TWO-LOOP SCALAR INTEGRALS 



Here, we report the identities that are used in the process of obtaining Eq. EH and 
Eq J191 they can be derived from the loop momentum translation invariant of the amplitudes. 
They are 
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f d D q 1 d D q 2 1 { q 2 - (g 2 - gi) 2 (gi • g 2 ) 2 - q\q 2 ^ _ 

J (2tt) d (2tt) d V \ D q\-m\ D{D - 1) J 
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with T> = ((g 2 — gi) 2 — m 2 )(q 2 — m 2 )(q 2 — m 2 ), and .D is the time-space dimension. In 
addition the two-loop vacuum integral 
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APPENDIX B: FORM FACTORS IN THE TWO-LOOP WILSON COEFFI- 
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